Incompressible Liquid, Stripes and Bubbles in rapidly rotating Bose atoms at v = 1 
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We numerically study the system of rapidly rotating Bose atoms at the filling factor (ratio of par- 
ticle number to vortex number) v = 1 with the dipolar interaction. A moderate dipolar interaction 
stabilizes the incompressible quantum liquid at v = 1. Further addition induces a collapse of it. The 
state after the collapse is a compressible state which has phases with stripes and bubbles. There are 
two types of bubbles with a different array. We also investigate models constructed from truncated 
interactions and the models with the three-body contact interaction. They also have phases with 
stripes and bubbles. 
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I. INTRODUCTION 

Recently rapidly rotating ultracold atoms in a trap 
have attracted considerable interest. Experimental stud- 
ies of moderately rotating atomic Bose gases have shown 
evidence for the formation of a vortex lattice^. Bose- 
Einstein condensation occurs even in the case without 
interaction but it is interaction which forms a vortex lat- 
tice. The number which parameterizes the population 
of vortices is given by v — N/Ny where N is the num- 
ber of bosons and Ny is the average number of vortices. 
For v > v* where v* ~ 6, the vortex lattice is stable^. 
For the high frequency regime v < i>* , numerical studies 
have predicted that the vortex lattice ought to melt and 
should be replaced by incompressible liquids which are 
closely related to fractional quantum Hall states^^. In 
that regime, v corresponds to the filling of the Landau 
levels in the rotating plane. At v = 1/2, the bosonic 
Laughlin stated is predicted^, while at v — k/2 (k is an 
integer > 2), states with a clustered structure emerge^. 
For k = 2 it is called the Pfaffian stated. Its fermionic 
counterpart is believed to be realized in the v = 5/2 
plateau^ and has been proposed as an excellent candi- 
date for a quantum computational device^. 

In typical Bose atomic gases, the interaction is short- 
ranged and can be modeled by a contact interaction 10 . 
Recently the condensation of Bose atoms with a strong 
magnetic dipole moment has been observed^. These 
atoms interact nonlocally through the dipolar interac- 
tion. The effect of the dipolar interaction on vortex lat- 
tices and incompressible liquids above has been studied 



for v = 1/2 1 



3/3M and v = 2±i. For these cases, 



a moderate dipolar interaction stabilizes the incompress- 
ible liquids. When the dipolar interaction gets strong, 
however, these states collapse. For the v — 1/2 case, 
the state after the collapse is identified as a stripe state 
in both Gross-Pitaevskii mean-field theory and micro- 
scopic theory 12 . As the dipolar interaction gets strong, 
the stripe state at v = 1/2 eventually shows a transition 
to a bubble stated 2 -. Similar transitions are also observed 
for v = 2^. A tendency of a spectrum to form stripes is 
mentioned also for v = 3/2^. These stripes and bubbles 



are closely related to those of two-dimensional fermions 
in higher Laudau level a 21 ' 22 i 23 . 

In this paper, we consider the finite system of Bose 
atoms at v = 1 (up to N = 12) which interact through 
the dipolar interaction and perform an exact diagonaliza- 
tion study to examine whether similar stabilization and 
collapse are induced. We observe the stabilization of the 
Pfaffian state similar to those of incompressible liquids 
at v = 1/2,3/2 by a moderate dipolar interaction. Fur- 
ther addition of the dipole moment induces a collapse 
of the incompressible liquid. To investigate the nature 
of the collapse, we take the periodic rectangular geom- 
etry. This geometry is suitable to address bulk prop- 
erties and transitions to other states. For the cases we 
have investigated up to N — 12, we observe four types of 
states: stripes, bubbles, intermediate states of two types 
of stripes, and intermediate states of stripes and bub- 
bles. These states are compressible and therefore sensi- 
tive to geometry. The periodic rectangular geometry is 
parameterized by an aspect ratio (< 1). We study the 
spectral flow of low-lying levels due to the change of the 
aspect ratio in detail. We also use the pair distribution 
function and the guiding-center static structure factor of 
ground states. They give a consistent picture of phases 
of the v = 1 Bose atoms. For the stripe states, the di- 
rection and the number of stripes depend on the aspect 
ratio and the interaction. For large aspect ratios, we 
observe that the increase of the strength of the dipolar 
interaction results in the formation of bubbles. Further 
addition of the dipolar interaction leads to a transition 
to another bubble state which has a different array. For 
small aspect ratios, stripe states seem to be favored even 
when the strength of the dipolar interaction is strong : 
we observe no evidence that a transition from stripes to 
bubbles occurs. 

To elucidate the nature of these transitions and real- 
ized states, we consider models constructed from trun- 
cated interactions. Any two-body interaction between 
particles in the lowest Landau level can be expanded in 
pseudopotential 15 V m ,m = 0, 1,2,3, each term of 
which corresponds to the projection to the two parti- 
cle states with definite relative angular momentum. For 
bosons, the odd terms vanish V m = 0, m = 1,3, 5, ■ ■ • . 
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For the dipolar interaction, all the even terms V m ,m = 
0, 2,4, • • • appear. We investigate the model only with 
the V% interaction ( "hollow-core" ) . It turns out that this 
model is always in a bubble phase for the aspect ratio 
above ~ 0.5. A certain amount of the two-body con- 
tact interaction Vq causes a transition to another bubble 
state with a different array. For small aspect ratios below 
~ 0.5, we observe the formation of stripes. 

These results are for two-body interactions. We also 
address the same issue to the three-body contact inter- 
action, for which the Pfaffian state is the exact ground 
stated. We investigate the effect of the Vi interaction on 
the model. Similar transitions to stripes and bubbles are 
observed also for this model. 

The organization of this paper is as follows. In Sec. [Til 
we give details of the formalism used to treat the bose 
problem on the periodic rectangular geometry. In Sec. 
IIII1 we show results for the dipolar model. In Sec lIVI is 
devoted to the V2 model. In Sec[V] we show results for 
the three-body contact interaction. Sec lVIl gives conclu- 



where g = 4:Trh 2 a s /m (a s is the s-wave scattering length). 
In typical Bose atoms, the interaction between particles 
is short-ranged and can be approximated excellently by 
Vq. Next, the electric or magnetic dipole interaction is 
given by 



\- Pi-Pj -3(n tJ j>i)(njj -pj) 

l<i<j<N 11 



Vdip — Cd 



where = (r, — Tj)/\ri — Tj\, the p-s are unit vec- 
tors which represent direction of the dipole moment. 
We assume that the dipole moments are parallel to the 
z-direction. The model with the dipolar interaction 
H = Vq + Vdip will be called as the dipolar model below. 
We expand H in pseudopotentiali5< V m ,m = 0, 1, 2, • • • . 
For bosons, only even m contribute. V m depends on the 
trap asymmetry l z /l. We assume that uj z is sufficiently 
large compared to lo so that we can take the limit that 
the thickness along the z-direction of the two-dimensional 
motion (in the xy plane) is negligible, l z /l — > 0. To first 
order in l z /l, the pseudopotential is given hy^ 



II. 



FORMALISM 



Vn = \ - 



2 h 2 a. 



il 2 l, 



The Hamiltonian describing N bosons of mass m in a 
rotating coordinate is 



H = 



N 1 N j 



i=l 

l<i<j<N 



(1) 



where the angular velocity is wz, the x-y trap frequency 
is ujq, and the axial trap frequency is lj z . The con- 
finement length in the rotating plane (the x-y plane) 
is I = ^/h/(mu>o) and the confinement length along 
the z-direction is l z = \J h/{muj z ). When co z is suf- 
ficiently large, the strong confinement along the z- 
direction freezes the motion along the z-direction in the 
ground state of the harmonic oscillator and the motion 
becomes quasi two-dimensional. For lj ~ uiq, the system 
is equivalent to the bosons of charge q in a magnetic field 
B = (2mui/q)z. When tu is sufficiently large and interac- 
tion is weak, we can treat the x-y motion to be restricted 
to the lowest Landau level. For N bosons spread over an 
area A, the characteristic parameter is v = N/Ny where 
Ny = (mLoA)/irh is the average number of vortices. In 
this paper, we focus our attention at v = 1. 

Let us next describe interactions between the atoms. 
The two-body contact interaction^ is given by 



Vo=9 £ 

l<i<j<N 



£ 3 (r 4 -t s ), 



(2) 



V, 



m>0 



7r(2m-3)!!C d 



m\2 r 



P 



(5) 



The strength of the dipolar interaction in the dipolar 
model V + Vdip is controlled by a = V 2 /V . We will also 
denote the potential with Vk ^ for k = i, VJ- = for 
k ^ i as Vi below. 

For realistic atoms, there are higher order terms of 
many-body interactions. The first of them is the three- 
body contact interaction 



^ b = c 3b 53 ^ 

l<i<j<k<N 



rj)8 3 (ri -r k ) 



(6) 



The exact zero-energy ground state of V^b is the Pfaffian 
state 16 : for the disk geometry in the symmetric gauge, 
it is given by (except for the gaussian factor of the lowest 
Landau level) 



pr 



Pf 



(7) 



K 3 



where Zi — Xj + iyi. This state is expected to have 
quasiparticles and quasiholes with nonabelian statistics 7 . 
The modular 5* matrices underlying the statistics have 
been studied in Refii 7 -. 

We consider several models with these interactions on 
the periodic rectangular geometry. We denote the sides 
of the rectangular as a and b. The translational invari- 
ance of the system gives a conserved pseudomomentum 18 
K = (K x ,Ky). The pseudomomentum runs over a Bril- 
louin zone containing N 2 points, where N is the great- 
est common division of N and Ny . K x , K y is measured 
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in units of 2nh/a and 2nh/b respectively. The states at 
(±K X , iK y ) are degenerate by symmetry so we may take 
only positive K x .K y . 

Next we recall the definition of some functions. The 
guiding-center static structure factor Sb(q) f° r a state 
|^>) is given by 



5 (q) =^(*|e iq ' (Ri - R ^|*) 



(8) 



hi 



where is the guiding center of the i-th particle(See 
ReiM for detail). Also, the pair distribution function in 
the real space G(r) is defined by 



G(r) 



N(N 



— (*|^«5(r-r 4 -r,)|*) (9) 



G(r) and So(q) will be used to study the nature of ground 
states. 
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FIG. 1: The squared overlaps of the two degenerate ground 
states and a nearly degenerate state for both the dipolar 
and the Vb — Vi model with the PfafRan states are shown 
for N = 10 at a/b = 1.0. The PfafRan states exist at 
K = (5, 5), (5, 0), (0, 5). The two degenerate ground states 
appear at K = (5,0), (0,5). The degeneracy is due to geo- 
metrical symmetry at a/b — 1.0. 



III. 



RESULTS FOR THE DIPOLAR MODEL 



In this section, we present results of exact diagonal- 
izations of the dipolar model Vb + Vdi P at v = 1. The 
strength of the dipolar interaction is controlled by a de- 
fined in the previous section. As discussed in Ref.— , the 
ground state at a = is an incompressible liquid and 
has a large overlap with the Pfaffian state. We add a 
moderate dipolar interaction (small a) and examine the 
stability of the incompressible liquid. 

We first fix the aspect ratio to be 1.0. For even N up to 
12, we observe that three nearly degenerate states with 
lowest energies appear at K where the Pfaffian states 
exist and are clearly separated from the rest of the spec- 
trum. 

In Fig. [H we present the squared overlaps of two de- 
generate ground states and a nearly degenerate state with 
the Pfaffian states for N = 10 as functions of a. At small 
a, the overlaps increase monotonically as a increases. 
The maximum of the overlaps is reached around a ~ 0.3. 
Thus a moderate dipolar interaction with a below 0.3 
stabilizes the Pfaffian states. Figure [2] shows the pair dis- 
tribution function G(r) of the Pfaffian state for N = 10 
and the ground state of the dipolar model for N = 12 at 
a = 0.3. Both states have very similar shapes of G(r). 
Except for regions near four corners on the rectangular 
boundary, G(r) is almost constant and isotropic. This 
property indicates incompressibility of the Pfaffian state. 
The G(r) of the ground state of the dipole model for 
N = 12 is slightly flatter than that of the Pfaffian state 
for N — 10. The nearly degenerate states in two other 
sectors of K have similar G(r). 

When the dipolar interaction gets strong, the overlaps 
show a very abrupt drop. It occurs around a ~ 0.4 for all 
the sectors as shown in Fig. [T] We observe similar stabi- 
lization and collapse occur in the region a/b = 0.3 ~ 1.0. 
The strength of the dipolar interaction a* where the col- 
lapse occurs vary with the aspect ratio. If we decrease 



G(r) 



2 r 




FIG. 2: Pair distribution function G(r) of various states at 
a/b = 1.0. (a) The Pfaffian state for N = 10 at K = (0,5). 
(b) The ground state of the dipolar model for N — 12 at 
a = 0.30, K = (6, 6). 



the aspect ratio, a* gradually increases and reaches to 
~ 0.5 at a/b = 0.3. We also confirm that similar stabi- 
lization and collapse occur for N = 6, 8 at a = 0.4 ~ 0.5 
for a/b = 0.3 ~ 1.0. 

We next investigate states after the collapse of the in- 
compressible liquid. In Fig. [31 we present the energy 
spectrum for N — 10 at a = 0.65 and a/b = 0.8. The 
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FIG. 3: Energy spectrum of the dipolar model for N = 10 at 
a = 0.65, a/b = 0.8 as a function of if = y/K£ + K%. The 
ground state at K = (0, 5) and states at K = (0, 3), (0, 1) are 
quasidegenerate . 



S (q) 




FIG. 4: Guiding-center static structure factor So(q) and pair 
distribution function G(r) of the ground state of the dipolar 
model for N = 10 at a = 0.65, a/b = 0.8, K = (0, 5). 



ground state at K = (0, 5) and the lowest energy states 
at K = (0,3), (0, 1) are quasidegenerate and clearly sep- 
arated from the rest of the spectrum. A characteristic 
wave vector that separates these quasidegenerate states 
is K* = (0,2). This suggests that these quasidegenerate 
ground states have a tendency to form a stripe state. 

Let Nd be the number of distinct quasidegenerate 
ground states, and let N s be the number of stripe. There 
are N 2 distinct values in the Brillouin zone, where N is 
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FIG. 5: Energy spectrum of the dipolar model for N = 12 at 
a — 0.6, a/b — 0.5. The ground state at K = (6, 0) and states 
at K = (6,3), (6,6) are quasidegenerate. The lowest energy 
states at K = (5, 1), (5, 2), (5, 4), (5, 5) form a quasidegenerate 
first-excited states. 



the greatest common divisor of N and Ny- If the trans- 
lational symmetry is broken in one direction, the length 
of the Brillouin zone of the one-dimensional superlattice 
is (2ttN)/(LNd) where L is the period of the superlat- 
tice. It must be 2ir/(L/N s ) where L/N s is the length 
per stripe. This yields N S N D = N. For v = l,N = N. 
The relation N s Nu = N yields the number of boson per 
stripe M = N/N s = N D . In the case above, N = 10 and 
Nd = 5, which gives N s — 2 and M = 5. 

To study properties of these quasidegenerate states, we 
calculate the guiding-center static structure factor So(q) 
and the pair distribution function G(r) of the ground 
state. We measure q x and q y in units of 2Tth/a and 2ith/b. 
In Fig. 21 we present S'o(q) of the ground state for N = 10 
at a = 0.65, a/b = 0.8, K = (0,5). There are peaks at 
±q* = (0,±2) and ±2q*. Other peaks at ±3q*,±4q* 
are due to the translational symmetry at v = 1. The 
peaks at ±q* = (0, ±2) indicate a strong density-density 
correlation in the ground state at this ordering vector. 
In Fig. 21 we also present G(r) of the state. The shape 
of G(r) is quite different from that of the Pfaffian state 
shown in Fig. [5] G(r) depends on the y-direction but 
has little dependence along the x-direction. Two peaks 
are seen along the y-direction in the unit cell, which is 
consistent with the peak of 5o(q) at q* = (0,2). So(q) 
and G(r) indicate that this ground state is a stripe state 
with two stripes lying parallel to the x-direction. This 
is consistent with the analysis of N s given in the previ- 
ous paragraph. So(q) an d G(r) of other quasidegenerate 
states at K = (0,3), (0,1) are similar to those of the 
ground state at K = (0, 5). 

In Fig. we present the energy spectrum for N = 12 
at a = 0.6, a/b = 0.5. The ground state at K = (6,0) 
and the lowest energy state at K = (6, 3), (6, 6) are 
quasidegenerate and clearly separated from the rest of 
the spectrum. These states are separated by a charac- 
teristic wave vector K* = (0, 3). The direction of stripes 
is parallel to the x-direction and the number of stripes 




FIG. 7: Energy spectrum of the dipolar model for N = 
10 at a = 1.0, a/b — 0.8. The quasidegenerate 
states are at K = (1,0),(0,1),(2,1),(1,2)(3,0),(0,3),(2,3),(3,2) 
,(4,1),(1,4),(4,3),(3 I 4),(5,0),(0,5),(4,5),(5,4),(2,5),(5,2). 



FIG. 6: So(q) and G(r) of the ground state of the dipolar 
model for N = 12 at a = 0.6, a/b = 0.5, K = (6, 0). 



is three. Fig. [6] shows 5o(q) and G(r) of the ground 
state. There are peaks at ±q* = (0, ±3) and at ±2q*. 
Other peaks at ±3q* are due to the translational sym- 
metry. We also present G(r) of the ground state in Fig. 
[5J There are three peaks along the y-direction. This is 
consistent with the peak of So (q) at q* . Quasidegenerate 
states at K = (6,3), (6,6) show similar So(q) and G(r). 
In Fig. we also observe a clear presence of nearly de- 
generate first-excited states formed by the lowest energy 
states at K = (5, 1), (5, 2), (5, 4), (5, 5). They are cre- 
ated from the quasidegenerate ground states by vectors 
ei = (1, 1), e2 = (—1, 1). A part of a less clear array of 
nearly degenerate second-excited states in Fig. [5] is also 
generated from the first-excited states by these vectors. 
Some other excited states in the second-excited array 
have the quantum numbers generated from (1, 0)(— 1, 0). 
The appearance of clearly separated low-energy bands is 
accounted for by low-energy particle-hole excitations in 
which a boson is moved from one stripe to another. These 
bands are almost flat, showing that such a hopping of the 
bosons between stripes is strongly suppressed. This band 
structure gives evidence of the stripe state. The band for 
excited states is not clear in Fig. [3] for N = 10. For 
N = 12, the stripe state with three stripes appears in 
the region a/b = 0.4 ~ 0.65. At a/b = 0.5, we con- 
firm that these states persist for a = 0.5 ~ 0.8. Simi- 
larly, we observe quasidegenerate ground states forming 
two stripes lying parallel to the y-direction in the region 
a/b = 0.7 - 1.0. For a beyond 0.6 at a/b = 0.7 - 1.0, 




FIG. 8: G(r) of the ground state of the dipolar model for 
N = 10 at a = 1.0, a/b = 0.8, K = (0, 5). 



quasidegenerate ground states persist to be present but 
their ordering is different from stripes (see below). 

Let us next turn to the region of larger dipole moment. 
In Fig. [71 we present the energy spectrum for N = 10 
at a = 1.0, a/b = 0.8. There are several quasidegener- 
ate states whose locations in the reciprocal space form a 
lattice. This means that these states have a tendency to 
form a bubble state. The primitive vectors are ei and 
e 2 . Let Nt be the number of bubbles per the unit cell. 
As derived in Refi22, the relation N^Nd — N 2 gives the 
number of bosons in a bubble M = N/Nj, = Nd/N. In 
this case, N = 10 and Nd — 50, which gives Nt, — 2 and 
M = 5. In Fig. [JJ one also sees quasidegenerate first- 
excited levels above the quasidegenerate ground states. 
A part of their quantum numbers are generated by vec- 
tors (±1, 0), (0, ±1) from those of the quasidegenerate 
ground states. The rest of them is the first-excited state 
at K where the quasidegenerate ground states appear. 
These almost flat bands are accounted for by low-energy 
particle-hole excitations in which a boson is moved from 
one bubble to another. In Fig. [51 we show the pair distri- 
bution function G(r) at K = (0, 5). There are two peaks 
at r = (0, 0) and r = (a/2, 6/2), which is consistent with 
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FIG. 9: Energy spectrum of the dipolar model for N = 10 at 
a = 3.0, a/b = 0.8. Every lowest energy state in the reciprocal 
space is quasidegenerate. 
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FIG. 10: G(r) of the ground state of the dipolar model for 
AT = 10 at a = 3.0, a/6 = 0.8, K = (0, 5). 



the number of A^. G(r) for other quasidegenerate states 
is similar. These give evidence that these quasidegen- 
erate states form a bubble state. The transition from 
the stripe state to the bubble state is observed around 
a ~ 0.8. For N = 8, bubble states are seen in the region 
a = 0.8~ 2.0, a/6 = 0.4 - 1.0. 

In Fig. [HI we present the energy spectrum for N = 10 
at a = 3.0, a/b = 0.8. The number of the quasidegen- 
erate ground states increases compared to the a = 1.0 
case. The lowest states for every K become quasidegen- 
erate. A clear presence of the band structure in Fig. [5] 
is readily interpreted by low-energy particle-hole excita- 
tions. These give evidence of a bubble state. In this 
case, N b = 1 and M = 10. In Fig. [TDJ G(r) of the 
lowest energy state at K = (0, 5) is presented. The den- 
sity concentrates near the corner of the unit cell, which 
is consistent with JV& = 1. Thus, the bubble state ob- 
served here is different from the bubble state for lower 
a. The transition between these bubble states occurs 
around a ~ 2.0 for N = 8, 10. 

As seen above, after the collapse of the incompress- 
ible liquid, ground states are not separated from excited 
states by a gap. This is characteristic of a compressible 
state. In general, its property is sensitive to geometry 



(the aspect ratio). Therefore we vary the aspect ratio 
for values of a and study the dependence of spectra on 
the aspect ratio. In Fig. QTJ we present energy spectra 
for N = 10 at a = 0.65, 1.0, 3.0. For a = 0.65, changes 
of the level structure accompanying level crossings are 
seen around a/b ~ 0.75,0.52,0.25. At a = 1.0, a struc- 
tural change is seen around a/b ~ 0.4. At a = 3.0, 
structural changes are seen around a/b ~ 0.7,0.55,0.3. 
Similar structural changes are observed for other N but 
the values of aspect ratios at which they occur depend on 
N. When structural changes are not sharp, intermediate 
states appear around these points as seen in Fig. QT] 

Let us illustrate a complicated behavior at a — 0.65 
in the region a/b = 0.5 ~ 0.7 where one sees many level 
crossings as in Fig. [TlTah As an example, we present 
the energy spectrum for N = 10 at a/b = 0.6 in Fig. [T2J 
The ground state appears at K = (5,5). The locations 
of quasidegenerate states in the reciprocal space do not 
form any lattice. For N — 12, similar states appear in 
the region of high aspect ratios a/b = 0.65 ~ 1.0 around 
a = 0.65. We interpret them as intermediate states be- 
tween stripe states and bubble states due to finite-size 
effects. The region for which such intermediate states 
appear should shrink as N gets large. 

Let us next turn to the highly anisotropic region (a/b 
below 0.5) of Fig. [TIT ah Energy spectrum and G(r) 
shows that stripe states appear in the region a/b = 
0.25 ~ 0.5. The quasidegenerate states in this region 
form a stripe state with three stripes lying parallel to the 
x-direction2^. 

Let us next illustrate some details of the level structure 
for a — 1.0 shown in Fig. fTTT b). The level structure of 
the region a/b — 0.2 ~ 0.4 is different from those for the 
bubble state seen in the region a/b = 0.4 ~ 1.0. The level 
structure of the bubble state in the region a/b = 0.4 ~ 1.0 
persists for a = 0.8 ~ 2.2. As seen in FigJTJJb), it disap- 
pears around a/b = 0.4 and changes to that of a stripe 
state. For other a's, this structural change occurs around 
a/b = 0.4 ~ 0.5. The direction of stripes is parallel to 
the ^-direction. 

In Fig. [UJc), we present energy spectra versus the 
aspect ratio for N — 10 at a — 3.0. It shows a change 
of the level structure around a/b ~ 0.7, 0.54. The bubble 
states disappear for lower aspect ratios. Stripe states 
appear at a/b below 0.54. 

In stripe states in Fig. HTT a).(b).(c) at the highly 
anisotropic region, the direction of stripes is always par- 
allel to the x-direction (the shorter axis). The number of 
stripes depends on N, a, and a/b. 

Let us next turn to the region of extremely low aspect 
ratios below 0.2. In this region, some strange behaviors 
are observed. For N — 8, 10, the overlap of quasidegen- 
erate ground states at a — with the Pfaffian states is 
large only for one sector of K . It drops rapidly when we 
move to a / but the other two sectors for quasidegen- 
erate ground states begin to have a large overlap with the 
Pfaffian states. The collapse of these two states occurs at 
much higher value of a than for nonextreme aspect ra- 
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FIG. 12: Energy spectrum of the dipolar model for N — 10 
at a — 0.65, a/b = 0.6. 
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FIG. 11: Energy levels versus the aspect ratio a/b for TV = 10. 
(a)a = 0.65, (b)a = 1.0, (c)a = 3.0. 



tios. In this region, the shapes of 5o(q) and G(r) for the 
Pfaffian states are quite similar to those of stripe states. 
They cannot be distinguished from a stripe state with 
four stripes (N = 8), five stripes (N = 10) lying parallel 
to the x-direction2£. 

Let us next present results for the isotropic case a/b = 
1.0 where geometric symmetry is enhanced. Also for 



G(r) 




FIG. 13: G(r) of the ground state of the dipolar model for 
N = 12, a = 0.6, a/b = 1.0 at K = (6, 6). 



this case, two stripe states appear for N — 10 in the 
region a — 0.4 ~ 0.8. States at K = (K\,K2) and 
K = (K2 , Ki ) are exactly degenerate due to geometri- 
cal symmetry present at a/b = 1.0. The ground states 
appear at K = (0,5) and K = (5,0). States at K = 
(1, 0), (3, 0), (5, 0), (0, 1), (0, 3), (0, 5) are quasidegenerate 
and separated from the rest of the spectrum. The states 
at K = (1, 0), (3, 0), (5, 0) form two stripes parallel to the 
y-direction while the states at K — (0, 1), (0, 3), (0, 5) 
form two stripes parallel to the ir-direction. Thus the 
ground state at a/b = 1.0 is a linear combination of 
these two stripe states. Similarly the ground state at 
K = (6, 6) for N — 12 is not a stripe state with a broken 
translational symmetry for one direction, nor a bubble 
state with a broken translational symmetry for two di- 
rections. In Fig. [T51 we present G(r) of the ground state 
at a = 0.6, K = (6,6). They indicate that the ground 
state has a strong density-density correlation along two 
stripes lying parallel to the x- and y-directions. We ob- 
serve that this shape of G(r) at a/b — 1.0 persists for 
N = 10, a = 0.4 - 0.8. 

Near the isotropic case, one expects that there may be 
a competence between stripe states with perpendicular 
directions. This is actually the case for N = 10 in a 
region around a/b ~ 0.8. The direction of the stripes 
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FIG. 14: G(r) of various ground states of the dipolar model, 
(a) N = 10, a = 0.49, K = (5, 0). (b) N = 10, a = 0.50, K = 
(0,5). 



changes as one varies a. There is a level crossing at 
a ~ 0.496 for a/b = 0.8. In the region a < 0.496, the 
ground state appears at K = (5,0). At a ~ 0.496, the 
lowest energy states at K = (5,0) and K = (0,5) are 
quasidegenerate. In the region a > 0.496, the ground 
state appears at K = (0,5). We present G(r) of the 
ground states for a — 0.49,0.50 in Fig. [14] (a),(b). In 
the region with a < 0.496, the quasidegenerate ground 
states at K = (5, 0), (3, 0), (1, 0) form a stripe state with 
two stripes lying parallel to the y-direction. At a = 0.50, 
the ground state at K = (0, 5) is an intermediate state 
between two stripes with perpendicular directions. At 
a = 0.65, the stripe state with two stripes lying parallel 
to the x-direction appears. 

Before ending this section, let us briefly discuss the 
V0-V2 model. Previously the collapse of the incompress- 
ible liquid for the two-body contact interaction Vo by 
the V2 interaction has been observed in the spherical 
geometr y 19 i 20 . In Fig. [H we present the squared overlaps 
of two degenerate ground states and a nearly degenerate 
state with the Pfafhan states for N = 10 at a/b — 1.0. As 
in the case of the dipolar interaction, a moderate V2 in- 
teraction stabilizes the Pfafhan state. It collapses around 
a ~ 0.25 as seen in Fig. [TJ We confirm that quasidegen- 
erate ground states with a tendency to form stripes ap- 
pear after the collapse. When we further increase a, two 
kinds of quasidegenerate ground states with characteris- 
tics to form bubbles appear as in the case of the dipolar 
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FIG. 15: Energy spectrum of the V2 model for N = 10 at 
a/b = 0.8. 




FIG. 16: G(r) of the ground state of the V 2 model for N = 12 
at a/b= 1.0, K = (6,6). 



model. 



IV. RESULTS FOR THE V 2 MODEL 

In this section, we present results for the V2 model. We 
present the energy spectrum for N = 10 at a/b = 0.8 in 
Fig. [I5l The lowest energy states at every point in the 
reciprocal space are quasidegenerate and separated from 
the rest of the spectrum. Also, in Fig. [T31 several flat 
excited bands are clearly seen. It is consistent with the 
presence of low-energy particle-hole excitations. These 
observations give evidence of a tendency to the formation 
of a bubble state. The pair distribution function G(r) of 
quasidegenerate states are similar to those of the bubble 
states of the dipolar mode at a — 3.0. Similar shapes of 
G(r) are observed for N = 8, 12 at a/b = 0.8 - 1.0. Fig. 
[Hlshows G(r) for N = 12, a/b = 1.0. 

We present energy spectra versus the aspect ratio for 
N = 10 in Fig. [T7J The spectral flow of the V 2 model is 
similar to that of the dipolar model at a above 2.0. As 
seen in Fig. [17l quasidegenerate ground states forming 
a bubble state persist for a/b = 0.5 ~ 1.0. The level 
structure changes to that of a stripe state around a/b = 







a 



0.2 0.4 0.6 0.8 1 

a/b 

FIG. 17: Energy levels versus the aspect ratio of the V2 model 
for N = 10. 



FIG. 19: The squared overlaps of the two degenerate ground 
states of the dipolar model with the ground states of the V2 
model are shown for N = 10 at a/b — 1.0. We present results 
for the two ground states at K = (5,0), (0,5). K = (0,5) 
and K = (5, 0) are exactly degenerate due to geometrical 
symmetry at a/b = 1.0. 




K 

FIG. 20: The squared overlaps of the two degenerate ground 
FIG. 18: Energy spectrum of the Vi model for N = 10 at states and a nearly degenerate state for the V3b + V2 model 
a/b — 0.35. The ground state at K = (0,5) and states at with the Pfafnan states are shown for N = 10 at a/b = 1.0. 
K = (0, 3), (0, 1) are quasidegenerate. The Pfafnan states exist at K = (5, 5), (5, 0), (0, 5). The 

two ground states of the V^b + V2 model appear at K = 
(5,0), (0,5). K = (0,5) and K = (5,0) are exactly degen- 
erate due to geometrical symmetry at a/b — 1.0. 



0.45. Fig. [T8l shows the energy spectrum at a/b = 0.35. 
The ground state at K = (0, 5) and the lowest energy 
states at K = (0,3), (0,1) are quasidegenerate. In Fig. 
[TBI clearly separated low-lying bands are seen, which is 
consistent with the existence of particle-hole excitations 
of a stripe state. The direction of stripes is always parallel 
to the x-direction. 

Fig[l9] shows the overlap of the ground state of the 
dipolar model with the ground state of the V2 model. 
The overlap is close to 1.0 for a beyond ~ 2.2. In the 
pseudopotential expansion, the dipolar model has higher 
order terms beyond V% but this result suggests that, for 
aspect ratios above ~ 0.5, quasidegenerate ground states 
for the dipolar model at large a coincide with those of 
the V2 model. 



V. RESULTS FOR THE THREE-BODY 
INTERACTION 

In this section, we present results for the Vsb + V2 
model. We set f3 = V2/C3&. We add a moderate V2 
interaction (small 0) and examine the stability of the 
Pfafnan states which are the exact zero-energy ground 
states at (3 = Q. In Fig. [501 we show the squared over- 
laps of nearly degenerate ground states at a/b = 1.0 with 
the Pfafnan states at each K where Pfafnan states exist. 
The overlaps decrease monotonically as (3 increases. The 
overlaps drop around /3~2.0xl0~ 3 . The collapse is not 
as sharp as that of the dipolar model. For other aspect 
ratios we have investigated, this collapse always occurs. 
The value of j3 where this collapse occurs increases as a/b 
decreases. 
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FIG. 21: Energy spectrum of the V 3b + Vi model for N = 10 
at (3= 1.0 x 10~ 2 ,a/& = 0.8. The ground state at K = (0,5) 
and states at K = (0, 1), (0, 3) are quasidegenerate. 



ratio, a stripe state reappear at a/b below 0.3. 

A bubble state appears at high aspect ratios for f3 = 
0.025 ~ 1.0. In Fig. [2H we present the energy spectrum 
for N = 10 at (3 = 0.1, a/b = 0.8. There are quaside- 
generate ground states arrayed in the lattice with the 
primitive vectors ei and e2, indicating a formation of 
bubbles. The presence of clearly separated bands in Fig. 
[22] gives further evidence of it. The G(r) of quasidegen- 
erate states is similar to that of the bubble state of the 
dipolar model presented in Fig. [5] At low aspect ratios 
a/b below 0.3, the quasidegenerate ground states remain 
to form a stripe state. The direction of stripes is always 
parallel to the rr-direction. The number of stripes varies 
for N and a/b. Intermediate states appear in the region 
a/b = 0.3 -0.5. 



VI. 



CONCLUSIONS 
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FIG. 22: Energy spectrum of the V-jb + Vi model for 
N = 10 at /3 — 0.1, a/b — 0.8. The quasidegenerate 
states are at K =(1,0),(0,1),(2,1),(1,2)(3,0),(0,3),(2,3),(3,2) 
,(4,1),(1,4),(4,3),(3,4),(5,0),(0,5),(4,5),(5,4),(2,5),(5,2). 



Let us next turn to states after the collapse. At high 
aspect ratios, we observe quasidegenerate ground states 
forming stripes. The stripe state changes to a bubble 
state around (3 ~ 2.5 x 10~ 2 . As we further increase 
/3, the ground state starts to have a large overlap with 
the ground state of the V% model at (3 above 1.0. The 
structure of the energy spectrum is similar to that of the 
V2 model. 

As an example of high aspect ratios, we present the 
energy spectrum for N = 10 at a/b = 0.8, f3 = 1.0 x 10~ 2 
in Fig. [21] The ground state at K = (0, 5) and states 
at K = (0,3), (0,1) are quasidegenerate. These states 
form a stripe state. The direction of stripes is parallel 
to the x-direction and the number of stripes is two while 
first-excited levels are not clearly separated in Fig. [5T] 
The pair distribution function G(r) for the ground state 
at K = (0, 5) is similar to that of the stripe state of the 
dipolar model presented in Fig. 0J When we decrease the 
aspect ratio, the stripe state changes to an intermediate 
state at a/b — 0.5. If we further decrease the aspect 



We have exhibited evidence that bosons with a dipole 
moment in the lowest Landau level at v = 1 have phases 
of incompressible liquid, stripes and bubbles based on 
exact diagonalizations up to N = 12. While a moderate 
amount of the dipolar interaction stabilizes the incom- 
pressible liquid, a further amount induces a collapse of 
it. The state after the collapse is a compressible state. 
Up to N = 12, it is a stripe state or an intermediate 
state between stripes and bubbles. The number and the 
direction of stripes are sensitive to the aspect ratio, the 
strength of the interaction and the number of particles. 
As the dipolar interaction gets strong, two kinds of bub- 
ble states appear for high aspect ratios. The transition 
between them is observed as the dipolar interaction gets 
strong. On the other hand, for small aspect ratios, stripe 
states remain to appear even when the dipolar interac- 
tion is rather strong. 

We have also considered the V2 model where the bosons 
do not have a hard-core. It has turned out that it does 
not have a phase of incompressible liquid. Its ground 
state is a bubble state for aspect ratios above 0.5. A cer- 
tain amount of the Vo interaction induces a transition to 
another bubble state a different array. We have observed 
that, for high aspect ratios, the quasidegenerate ground 
states for the dipolar interaction at the strong limit have 
large overlaps with those of the V2 model. For low aspect 
ratios below 0.5, quasidegenerate ground states of the V2 
model are identified to be a stripe state. 

We have also studied the model where the bosons inter- 
act through the three-body contact interaction for which 
the Pfaffian state is the exact ground state. It shows sim- 
ilar properties as those of the dipolar model : our results 
give evidence that the model has phases of stripes and 
bubbles. 

Finally we remark that the geometry treated in this 
paper is a subset of the general geometry of torus where 
the period vectors are not necessarily perpendicular. 
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Note Added. While we were completing this work, 



a preprint which has a partial overlap with this paper 
appeared^. Energy levels which show degenerate ground 
states which are characteristic of the Pfaffian and stripes 
are observed at a strength of the dipolar interaction. 
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